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The derivation of the macroscopic transport 
equations of hydrodynamics and the computa- 
tion of the associated coefficients for systems de- 
scribed at the microscopic level by Hamilton's 
equations of classical mechanics is a central prob- 
lem of non-equilibrium statistical physics. The 
periodic Lorentz gas provides an example where 
this program can be achieved and Pick's law of 
diffusion established^ Furthermore, by tweak- 
ing the system's geometry so that tracer particles 
hop from cell to cell at nearly vanishing rates, 
memory effects disappear and the dynamics of 
tracers is well approximated by a continuous time 
random walk^. In this regime, the diffusion coef- 
ficient takes on a simple limiting value, given by 
a dimensional formula^, by which we mean that 
its expression reduces to the square of the length 
scale of the cell separation multiplied by the hop- 
ping rate. Similarly, it was found that heat trans- 
port in systems of confined hard disks with rare 
interactions reduces to a stochastic process of 
energy exchanges which obeys Fourier's law of 
heat conduction, with the coefficient of heat con- 
ductivity also given by a dimensional formula'' 
where the timescale is that of energy exchanges 
among particles in neighboring cells. Here we ex- 
tend these findings to mechanical systems of con- 
fined hard spheres, whose stochastic limit was al- 
ready studied elsewhere^, and review in some de- 
tails the analogy with the problem of mass trans- 
port in the periodic Lorentz gas. 



I. INTRODUCTION 

A long-standing problem in non-equilibrium statisti- 
cal mechanics has been to provide a derivation from 
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first principles of Fourier's law of heat conduction in 
insulating materials in the framework of Hamiltonian 
mechanicJiSl. A centerpiece of the puzzle is to identify 
the conditions under which scale separation occurs so the 
time evolution at the microscopic level can somehow be 
reduced to a hydrodynamic equation at the macroscopic 
level, which is an especially challenging problem for in- 
teracting particle systems. 

As emphasized in earlier papers^" ^, the following two- 
step program, which consists of (i) identifying an in- 
termediate level of description — a mesoscopic scale — 
where the Newtonian dynamics can be consistently ap- 
proximated by a set of stochastic equations, and (ii) 
subsequently analyzing the statistical properties of this 
stochastic system and computing its transport properties 
in the hydrodynamic scaling limit, can be successfully 
achieved in a class of chaotic billiard systems composed 
of many hard disks trapped in a semi-porous material 
which prevents mass transfer and yet allows energy trans- 
fer through elastic collisions among neighboring disks. 

Under the assumption that collisions among moving 
disks are rare compared to wall collision events, the global 
multi-particle probability distribution of the system typ- 
ically reaches local equilibrium at the kinetic energy of 
each individual particle before energy exchanges proceed. 
This mechanism naturally yields a stochastic description 
for the process of energy exchanges in the system. In a 
subsequent papei^, it was shown that the same reduction 
applies to systems of confined hard spheres. 

Our purpose in this paper is twofold. Our main objec- 
tive is to establish a comparison between the process of 
energy transfer in models of interacting hard spheres with 
local confinement rules, on the one hand, and the diffu- 
sive motion of tracer particles in low-dimensional billiard 
tables such as the finite-horizon periodic Lorentz gas, on 
the other. Specifically, we show that the two systems, un- 
der equivalent assumptions of local equilibration, whose 
accuracy can be precisely controlled by tuning the sys- 
tems' parameters down to a critical geometry, are both 
amenable to stochastic descriptions in the form of mas- 
ter equations, whether for the distribution of energies in 
the case of the former systems, or that of mass in the 
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latter. Taking the hydrodynamic limit of these stochas- 
tic systems, we obtain explicit values of the transport 
coefficients of heat conduction and diffusion respectively, 
which turn out to share the remarkable property that 
they are given by simple dimensional formulae, i. e. by 
the square of the mesoscopic length scale multiplied re- 
spectively by the rates of energy or mass transfer. 

Turning back to the billiard dynamics, we address our 
second objective, which is to show that, under the as- 
sumption of separation of two characteristic timescales, 
one associated with the local dynamics, the other with 
energy transfers, the process of heat transport in three- 
dimensional billiard systems of confined hard spheres is 
well-approximated by the corresponding stochastic pro- 
cess of energy exchanges. We do so by considering the 
Helfand moment of thermal conductivity^^ and compute 
the linear divergence in time of its mean squared change. 
Plotting our results as functions of the system's sizes for 
different parameter values, we compute the infinite-size 
extrapolations to obtain the heat conductivity and ob- 
serve a very good convergence to the value obtained for 
the stochastic systems for parameter values where an ef- 
fective separation of timescales is observed. 

The paper is organized as follows. The problem of mass 
transport in a spatially periodic billiard table is consid- 
ered in Sec. |TTj Starting from the pseudo-Liouville equa- 
tion for the billiard dynamics, we derive a continuous- 
time random walk which describes the stochastic jumps 
of particles across a lattice, and analyze its transport 
properties, comparing it to that of the billiard. The 
same procedure is applied in Sec. |III| to two-dimensional 
billiards of confined hard disks. In Sec. IV we turn to 
a three-dimensional billiard system and discuss our nu- 
merical results. Conclusions are drawn in Sec. IVl 



II. A CONTINUOUS-TIME RANDOM WALK 
APPROACH TO MASS TRANSPORT 

We consider for the sake of the example a periodic ar- 
ray of two-dimensional semi-dispersing Sinai billiard ta- 
bles in the form of square billiard cells bounded by flat 
walls of sizes I, and connected to each other by small 
openings of relative widths S. All the cells are identical 
and contain circular obstacles arranged so that there are 
no periodic trajectory that avoid these obstacles. Inde- 
pendent pointwise tracer particles move across this array, 
performing elastic collisions on the obstacles and walls. 
An example is shown in figure [T] The diffusive properties 
of this model were studied in some details in Ref. [T^ . 
There the attention was on memory effects and their im- 
pact on the value of the diffusion coefficient. Here we 
will focus on the derivation of the kinetic prediction of 
this coefficient, which yields the aforementioned dimen- 
sional formula, disregarding the memory effects. The 
conclusions are similar to those obtained by Machta and 
Zwanzi^ in the framework of the periodic Lorentz gas. 
The continuous-time approach we present below is some- 




FIG. 1. Periodic billiard table on a square lattice with a 
typical trajectory. Here the central disks in the initial and 
final cells are color-filled. One considers the process of mass 
transport across the periodic cells. The model has several 
parameters in terms of which a dimensional formula of the 
diffusion coefficient is computed. 



what similar to that of Zwanzi^. 

The phase-space configuration of this system refers to 
a single tracer particle and is specified by the triplet 
{n, r,v}, where n = (n^jUy) G denotes the lattice 
index of the cell where the tracer is located, r, its posi- 
tion within the cell and v its velocity. 

Let p{n, r, v, t) denote the probability distribution of 
the tracer. Its time evolution is determined by the 
pseudo-Liouville operator^, which comprises four differ- 
ent types of contributions, corresponding to as many dif- 
ferent types of events: 

1. Free advection of the tracer inside the cell is ac- 
counted for by the term — v • d^', 

2. Collisions of the tracer with any of the circular ob- 
stacles inside the cell are determined by the opera- 
tors K^'^^ , where d refers to a specific disk of radius 



Pd at position q' 



d . 



if('^)p(n,r,v,t) [ de(e-v) 

Jev>0 

X 6{r — — pe)p{n, r, v — 2e(e • v), i) 

-(5(r-q'* + /?e)p(n,r,v,i) ; (I) 



Collisions of the tracer with any of the flat walls 
of cell n, with operator W^^\ where j takes on 
the values j = 1, . . . , 4, corresponding to the right, 
bottom, left and top walls: letting and Vx denote 
the position and velocity components of the tracer 
along the horizontal axis, wall #1 is the right wall 
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where A is the area of the bilhard, e. g. A = l'^—'n'{pi+P2) 
in the case of figure [T] (where pi denotes the radius of the 
disk at the center of the cell and p2 that of the disks at 
the cell corners). 

Plugging Eq. ([t]) into ([6|, the pseudo-Liouville equa- 
tion Q reduces to the following continuous-time random 
walk 

with similar expressions for W^'^\ W'-^\ and W'-^'^; dtP'-^'"^\n,t) — ^ ^^[P^'°'^^(n + e^, — P^''''^^(n, t)] , 



at position 1/2 so that, if the wall is continuous with 
respect to the vertical axis, i. e. (5 = 0, the collision 
operator acts according to 

M^(iV(n,r,v,i) = |«,|5(r,~ V2) (2) 
'-Vx)p{n, r, -Vx, Vy, t) - e{vx)p{n, r, v, t) , 



4. Jump events between neighboring cells take place 
when a wall collision event occurs at a position 
where the wall is actually open, which amounts 
to subtracting the action of the wall operator ([2| 
above, keeping track of the cell indices: using the 
same set of four indices as above for every wall and 
assuming that wall #1 has a slit of width 6 centered 
about y — 0, the jump operator jj^' is 

J^^V(n, r, V, t) = IvMr. " l/2)esiry)e{-v^) (3) 
X p{n -h (1, 0), r, V, t) - p(n, r, -v^,Vy, t) 

where Qs{x) = 1 if < S/2 and otherwise. 

Collecting these terms together, we obtain the time evo- 
lution of p(n, r, V, t), given by the pseudo-Liouville equa- 
tion: 

dtp{n,r,v,t) = (4) 

{-V • 9r + Ed ^^'^ + } p(n, r, v, t) 

+ E,4'V(n,r,v,i). 

The terms on the right-hand side of this equation are 
grouped so as to distinguish the local terms, which do 
not act on the cell index in the distribution, from the 
non-local ones, i. e. the jump terms, which act on the 
cell index. Also note that the velocity amplitude is in- 
variant under every term of this equation and is therefore 
a conserved quantity which plays no role once it has been 
fixed by the initial condition. 

It it an immediate consequence of the crgodicity of the 
local dynamics that the following distribution. 



P(''=i)(n,t) = J dr J dvS{v^ ~vl)p{n,r,v,t) 



(5) 



is invariant under the local terms of the pseudo-Liouville 
equation We thus refer to P('''i)(n,i) as the lo- 

cal equilibrium distribution. Its time evolution occurs 
through jump events only : 

5^p(lcq)(n ^ 

E, / dr / dvS{v^ - v'o)4'^p{n, r, v, t). (6) 

In order to close this equation for p('°i-' (n, t), we make 
the local equilibrium approximation, 

p{n, r, V, t)^-^6{v'- v^P^"'^^ (n, t) , (7) 



J 

(8) 

where e {(±1, 0), (0, ±1)}. 

In the continuum scaling limit, f — In, ^ — > 0, wq ^ 
this reduces to a diffusion equation. 



dtPif,t)^VV^P{f,t), 



with diffusion coefficient 



2? = ;2^w, 



(9) 



(10) 



here written in terms of the frequency i^^"'-' of a jump 
event in a specific direction, identical for all four direc- 
tions. 



ttA' 



(11) 



Equation ( 10 ) is what we refer to as a dimensional 
formula: the diffusion coefficient associated with the 
continuous-time random walk process ([8| is the product 
between the scale of displacements squared and hopping 
rate. The fact that it is an exact property of the stochas- 
tic process is remarkable. 

As far as the transport properties of the billiard table 
go, Eq. ( 10 1 is a mere approximation: its validity relies 
on the local equilibrium approximation Q. In contrast. 



Eq. (Ill is an exact result for the billiard dynamics as 
well — valid for all parameter values — as it relies solely on 
the ergodicity of the billiard tabl^^. 

The closure approximation ([t]) is expected to be exact 
in the limit of vanishing window widths, (5 — > 0, whereby 
the hopping rate (111 diverges and the diffusion coeffi- 
vanishes. 



cient ([10| vanishes. In practice however, we may expect 
p{n, r, V, t) to converge to the local equilibrium distribu- 
tion p('°i) (n, t) between successive jumps so long as the 
typical number of disk collision events within a cell is 
large. This can be verified numerically, as shown on the 
top panel of figure [2j These numerical results were al- 
ready reported in Ref. [T2]. Similar results are equally 
obtainable for three dimensional billiard systems such as 
the three-dimensional periodic Lorentz gas^. 



III. HEAT TRANSPORT IN TWO-DIMENSIONAL 
CONFINING BILLIARDS 

Similar considerations apply to models of heat trans- 
port with mass confinement^. Classes of such mod- 
els were initially introduced by Bunimovich et. al in 
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FIG. 2. (Top) Diffusion coefRcient as measured from tiie 
mean-squared displacement divided by the number of jumps 
for the Sinai billiard table shown in figure [l] (Bottom) Corre- 
sponding decay time. The dots show numerical measurements 
and the red line is the formula 




FIG. 3. Example of a two-dimensional billiard table consist- 
ing of many moving hard disks (colored), each trapped within 
confining walls (black disks) , here arranged in a square tiling. 
The centers of the moving disks are bound to the areas delim- 
ited by the exterior intersections of the black circles around 
the fixed disks (the solid broken lines sample their trajecto- 
ries). The parameters (radii of the black circles p and moving 
disks pm) are so chosen that (i) black circles overlap, so the 
moving disks are confined, and (ii) collisions are allowed (such 
as in the upper right corner), whereby energy exchanges take 
place. The mobile disks are color coded from blue to red with 
growing kinetic energies. 



Ref. [16 . There, the authors proved the ergodicity of 
two-dimensional billiard tables consisting of an arbitrary 
large number of unit cells placed side by side, each con- 
taining a single disk trapped within the cell's boundaries, 
but in such a way that collisions may still take place be- 
tween disks belonging to neighboring cells. 

Figure [s] shows an example of such a systenP^. In 
Ref. [7], we introduce the following parametrization of 
the model in terms of two parameters: p, which char- 
acterizes the timescale of wall collision events, and pni, 
which characterizes that of binary collisions. 

On the one hand, the dynamics within an isolated cell 
of width I boils down to the motion of a point particle 
in an area bounded by the exterior intersection of four 
disks of radius p, Z/2 < p < 1/^/2. In the absence of 
interaction with neighboring particles, the mean free path 
of a particle is given by ^ = 7r|Sp|/|9Sp|, where Bp and 
dBp denote respectively the area and perimeter of the 
billiard cell. The corresponding wall collision timescale 
is obtained by multiplying the mean free path by the 
speed of the particle. 

Interactions among neighboring particles on the other 
hand take place provided the radius of the moving par- 
ti cles Pm is larger than a critical value, Pm > Pc = 
\/ — P l'^. The corresponding timescale can be com- 



puted and shown to diverge with (p^ — Pc) Letting 
i/ON){T) and v^^^T) respectively denote the frequencies 
of wall collision and binary collision events measured at 
equilibrium temperature T, the following separation of 
timescales is assumed, 



^(B)(r)«^(W)(y)^ 



(12) 



which occurs when — )■ p^ . 

Under this assumption, we showed^ that the heat 
conductivity of spatially extended billiard systems, which 
are infinite size limits of systems such as depicted in 
Fig. [sj reduces, up to a dimensional factor P , to the 
frequency of binary collisions. 



«;(B)(T) = ?2j,(B)(y), 



(13) 



Furthermore the heat conductivity scales with the ther- 
mal speed, k(^) (T) - ^/T. 



Although Eqs. ( 10 1 and (^13| are very much alike, the 



derivation of the latter is much more involved than that 
of the former. 

Proceeding as in Sec. |ll] the first part of the pro- 
gram, which is to reduce the pseudo-Liouville equation 
describing the time evolution of the billiard system to a 
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continuous-time stochastic process of energy exchanges, 
follows closely along lines which led from Eqs. Q to 
Q . Indeed, considering the time evolution acting on the 
A^-particle phase-space distribution p7v({ri, v^}, t), one 
readily notices that the local equilibrium distribution 



N 



(14) 



J Y[ dradvQ(5(ea - mvl/2)pN{{r,,v,},t) 



is left unchanged by the advection and wall collision 



terms. Therefore, and provided binary collision events 
are rare with respect to wall collision events, we can make 
a closure approximation similar to Eq. ([7]) to obtain the 
time evolution of the local equilibrium distribution in the 
form of a stochastic process. 



The result is a master equation which describes the 
time evolution of the N cell system with energy variables 
{ei, . . . , eAf} in terms of energy exchanges between neigh- 
boring cells at respective energies and ei, of amount rj, 
specified by a stochastic kernel W: 



N 



a,b=l 



dtP^^'''^'{ei,...,eN,t) ^ / V7(ea 77, £6 - 77|ea, eb)-Pi) (• ••,£« + ?7,...,efc-?7,...,0 



W{ea, eb\ea- ??, eb + ri)pj^'''^\. .., e^, ...,£&,..., , 



(15) 



where the expression of W can be obtained by direct com- 
putation of the collision integrals, which, after rescaling 
the time variable to the units of the frequency of binary 
collision events and thus absorbing all the parameters 
into the time scale, yields the universal functiorP: 



W(ea,eb\ea-r],<^b + r/) 



K 



X < 



(^) 
( 
( 



ib 



K 



(16) 

whose defition intervals correspond respectively to — 6^ < 
T] < — max(ef, — Ca, 0), — max(e{, — Ca, 0) < 77 < 0, < 77 < 
max(ea — eb, 0), and max(ea — ef,, 0) < 77 < ea. Here KOm) 
denotes the complete elliptic intergal of the first kincpS. 
A system of TV isolated cells whose time evolution is 



specified by the master equation ( 15 ) reaches a micro- 



canonical equilibrium state whose total energy can be 
parametrized in terms of the temperature according to 
£1 -f . . . + e^v = NT. The corresponding energy exchange 
frequency is 



(17) 



un{T) = VT 1 + 0{1/N) , 



whose infinite-size limit is simply i^iT) = 
limi\i^rxi vn{T) =^ \/T (in the chosen units of time). 

An expression of the heat conductivity of the system 
described by Eq. (151 is obtained by considering the 



Einstein-type relation satisfied by the variance of the as- 
sociated Helfand moment, which measures the spread of 
energy as a function of time. 

For the system of N energy cells aligned along a one- 
dimensional ring, the Helfand moment is defined accord- 



ing to 



N 



(18) 



where ei(t) is the state of the energy at site i at time t. 
This quantity evolves in time by discrete steps, when en- 
ergy exchanges occur. Let {r„}„gN denote the sequence 
of times at which successive energy exchanges take place. 
Assuming cells i and i -I- 1 exchange some amount of en- 
ergy at time t„, we can write the corresponding change 
in the Helfand moment as ei(T„ — 0) — ei(r„ -|- 0). 

Computing the mean squared change in the Helfand 
moment as a function of time, we obtain an expression 
of the thermal conductivity according to 



k(T) = lim Kjv(r), 

N^oo 



(19) 



where we defined the finite N conductivity to be 
kn{T) 



— ^ lim ^([i/^(r„)-i/iv(ro)]2 

(20) 

In Ref. [8 , it was argued that the infinite N limit of 
this quantity is determined by the static correlations only, 
yielding the result 



KN{T) = fVT\l + 0{l/N) 



(21) 



Comparing Eqs. (17 1 and (21), we obtain the an- 
nounced result 

k{T) - 1\{T), (22) 
which is an exact result for the stochastic system evolved 



by the master equation ( 15 1 and does not make explicit 
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use of the form (16), except for some symmetric^. The 



corresponding resuh ( 13 1 for the bilhard dynamics is ob- 



tained by plugging back the proper timescale of binary 
collisions and letting — > Pc so that the separation of 
timescales ( 12 ) is effective. 



Unlike mass transport in Sinai billiard tables for which 
the transport equation ^ follows directly from the 
continuous-time random walk so that there only re- 
mains the problem of comparing the diffusion coefficient 
of the billiard to the dimension formula ( 10 1 in the ap- 



propriate parameter regime, the problem of computing 
the heat conductivity associated with heat transport in 
billiard systems of many confined particles proceeds in 
two separate steps. 

Having carried out the reduction of the billiard's 
pseudo-Liouville equation to a master equation for a 
stochastic energy exchange system, the first step is to es- 
tablish the cancellation of dynamical correlations in such 
a system, which yields the identity limTv-foo i^n I^n = 
This is however a delicate result and requires in-depth 
knowledge of the spectral properties of the master equa- 
tion (15). The approach we took in Ref. .8^, though 
equivalent, uses different techniques, based on analyz- 
ing the first few orders of the gradient expansion of the 
kinetic equation to obtain the expression of the heat cur- 
rent in terms of the local temperature gradient in a non- 
equilibrium stationary state, i. e. Fourier's law. 

The second step of the program is to go back to the 
billiard dynamics and investigate the convergence of the 
heat conductivity to the binary collision frequency in the 
limit pra -> Pc- 

For two-dimensional billiard systems such as the ones 
shown in Fig. [Sj the agreement was found to be rather 
satisfactory^. For the three-dimensional billiard systems 
we turn to now, this agreement is even better. 





IV. THREE-DIMENSIONAL BILLIARD FOR HEAT 
TRANSPORT 

Consider a three-dimensional billiard of confined hard 
spheres such as shown in Fig. [4] The reduction of the 
pseudo-Liouville equation governing the phase-space evo- 
lution of probability densities p7v({ri, v^}, t) to a master 
equation similar to Eq. (15 1 was already carried out in 
Ref. [g. 

The corresponding kernel, in the appropriate time 
units, is found to have the universal form 



1 



(23) 



whose definition intervals correspond respectively to 
~eb < ij < -max(eh - 6^,0), -max(eb - 6^,0) < < 
max(ea - 0), and max(ea - ef,, 0) < 77 < ea. 



FIG. 4. (Top) Hard-sphere particle trapped in a cuboid cell 
with cylindrical edges. (Bottom) A system made out of many 
copies of such cells which form a spatially periodic structure. 
The cells are semi-porous in the sense that particles are pre- 
vented from escaping, and can yet partially penetrate into 
the neighboring cells, thus allowing energy transfer through 
collisions among neighboring particles. As with the two- 
dimensional case, the likelihood of binary collision events can 
be controlled by tuning the geometry of the cell. 



Extensive numerical investigations of the master equa- 
tion ( [l5| with the stochastic kernel (23) were presented 
in Ref. I9j , supporting with high precision the validity of 
Eq. (22). Here we focus on the billiard dynamics and 



consider the mean squared change in time of the Helfand 
moment associated with the distribution of energy in bil- 
liard systems formed by one-dimcnsional lattices of bil- 
liard cells such as shown in Fig. |4] 

Let {Tn}nez denote the times at successive binary col- 
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lision events. As we now need to account for the motion 
of particles within their respective cells, there are two 
types of contributions to changes in the Helfand moment 
between two successive binary collision events. The lead- 
ing contribution arises from the energy exchanges which 
take place at binary collision events. Thus, when a binary 
collision occurs between particles j and fc, the Helfand 
moment changes by the amount [a;j(T„) — Xfc(r„)][ej(T„-|- 
0) — ej(T„ — 0)], where a;j(r„) and a;fe(r„) denote the po- 
sitions of particles j and k along the direction of the 
lattice spatial extension at time r„. The other contri- 
bution to the Helfand moment arises from the advec- 
tion of particles within their respective cells according to 

Y^ai^ai'^n) - (t„_ i )] (t„_ 1 ) . 

The range of allowed parameter values is identical to 
the two-dimensional case: 1/2 < p < l/V^, and pc < 
Pm < P, where p^ = p^ - P/A. 

Taking Z = 1, we fix p = 0.50 so pc = and use the six 
parameter values pm = 0.20, 0.25, . . . , 0.45. For each one 
of them, we fix the total energy to be S = 37V/2 (T = 1) 
and vary the system size from A'^ = 3 and up to TV = 100 
cells. We typically run 10'^ trajectories, each for a dura- 
tion of 10'^ X N units of the average time between collision 
events. We compute: (i) the Helfand moment versus time 
and take the mean and standard deviation of the trans- 
position of Eq. ( 20 ) to obtain the conductivities k]y '^ (ii) 

,(W) 
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FIG. 6. Ratio between the heat conductivity of the infinite 
length system and the corresponding collision frequency be- 
tween neighboring particles as a function of the parameter pm 
(p = 0.50 and pc = 0). The error bars show the widths of the 
computed 0.95 confidence intervals. 



sional formula (13 1 for values of pm as large as 0.30. In 
Fig. [TI the details of the fitting procedure used to obtain 
the values of Kqo and z^oo are shown for the different 
parameter values as functions of N . 



the binary and wall collision frequencies v)^ and 

by directly averaging their numbers with respect to time. ^- CONCLUDING REMARKS 



We then use linear fits in 1/iV for > 10 of k)^ and 

v^^^ and 



v^^^ to obtain the extrapolation v^r^^ 
= limjv^or, K^^^ 



lim 



N- 



N 



Here and below, explicit tempera- 



ture dependencies are dropped. 
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FIG. 5. Ratio between the binary and wall collision frequen- 
cies as a function of the parameter pm for N — 25 {p = 0.50, 
pc = 0). 



Measurements of the ratio between and z/]y , by 
which we can assess the effectiveness of the separation 
of timescales (12), are displayed in Fig. [s] for A^ = 25 
(other system sizes yield similar values). The results of 

(a) (a) 

the computations of Koo /voo for the different parame- 
ter values Pm are displayed in Fig. [6] and are found to be 
in very good agreement (within two digits) of the dimen- 



By controlling the rate at which a tracer hops from 
cell to cell in a Sinai billiard table or that of interaction 
among neighboring disks or spheres in a high-dimensional 
billiard with local confinement rules, one identifies a limit 
of vanishing rate where the complicated phase-space dy- 
namics are replaced by stochastic processes which ac- 
count for the transport of mass in the first case, or (ki- 
netic) energy in the second. Furthermore, the trans- 
port coefficients of these stochastic processes have the 
same simple dimensional expression, given by the length 
scale of transfers squared multiplied by the corresponding 
rates. 

As emphasized in this paper, the dimensional formulae 
we obtained for the transport coefficients of the models 
of mass and heat transport we have considered are but 
two faces of the same coin. Indeed, in both cases the ac- 
curacy of our approximation of the transport coefficients 
of the billiards in terms of hopping or collision rates relies 
on the efficient separation of two timescales, namely the 
timescale of local collision events must be much shorter 
than that characterizing the transfer of mass or energy. 
In other words, the relaxation to local equilibrium pre- 
cedes mass or energy transfers. 

We can in fact view the notion of relaxation to local 
equilibrium as a low-dimensional transposition of that of 
local thermal equilibrium, which is at the heart of many 
theories involving hydrodynamic scaling limits and typi- 
cally assumes a large number of degrees of freedom^*^. In 
our billiards, the relaxation to local equilibrium occurs 
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FIG. 7. Heat conductivities kJ^' and binary collision frequencies v)!^' measured as functions of the system size for different 



(B) 



values of the parameter pm (p = 0.50). The sohd curves correspond to the linear fits of and v'^^' as functions of 1/N . The 



intercepts yield the infinite size estimates /t^' and u[ 



,{B) 



on the constant energy surface of a single particle. But 
one could instead consider systems of trapped gases, in 
which case the relaxation to local equilibrium would in- 
volve transfers of energy among the particles in the same 
trap. In such a case, provided relaxation to local equi- 
librium takes place on timescales much shorter than that 
of energy transfer between neighboring traps, a similar 
dimensional expression of the heat conductivity of the 
corresponding stochastic process in terms of the prod- 
uct of length scale squared and rate of energy transfer 
would yield an accurate approximation of the transport 
coefficient of the billiard system. 

We end with a remark concerning the closure ([7| which 
relies on the ergodicity of the local dynamics on the sur- 



face of constant energy. As noted already by Zwanzi^, 
Sinai billiards cannot be replaced by polygonal ones, for 
which the notion of ergodicity is weaker since the ve- 
locity directions take on values in a discrete set. In 
energy exchange processes, however, ergodicity of the 
many-particle billiard may be restored through the sole 
interaction among neighboring particles. As shown in 
Ref. [20] . an elastic string- type interaction between par- 
ticles trapped in polygonal boxes provides a simple model 
of a system which, on the one hand, cannot be accurately 
described by a master equation similar to Eq. (15), even 



when interactions are rare, but on the other hand, has 
a well-defined heat conductivity which is well approxi- 
mated by a dimensional formula. An understanding of 
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the transport properties of this system beyond the Boltz- 
mann hypothesis remains to be elucidated. 
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